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"CHAIN POOLING" MODEL SELECTION AS DEVELOPED 
FOR THE STATISTICAL ANALYSIS OF A ROTOR 
BURST PROTECTION EXPERIMENT 

by Arthur G, Holms 
Lewis Research Center 

SUMMARY 

A statistical decision procedure called chain pooling had been devel- 
oped for model selection in fitting the results of a two-level fixed effects 
full or fractional factorial experiment not having replication. The basic 
strategy included the use of one nommal level of significance for a pre- 
limmary test and a second nominal level of significance for the fmal test 

The subject has been reexamined from the point of view of using as 
many as three successive statistical model deletion procedures m fitting 
the results of a single experiment. The new investigation consisted of 
random number studies mtended to simulate the results of a proposed 
aircraft turbine engme rotor burst protection experiment. As a conser- 
vative approach, population model coefficients were chosen to represent 
a saturated 2 ^ eiqjeriment with a distribution of parameter values un- 
favorable to the decision procedures. 

Three model selection strategies were developed, namely, (1) a 
strategy to be used when the e3q)ermienter anticipates a large error 
variance (coefficients of variation in the neighborhood of 66 percent), 

(2) a strategy to be used when the experimenter anticipates a small 
error variance (coefficients of variation of 4 percent or less), and (3) 
a security regret strategy to be used m the absence of such prior 
knowledge. 
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INTRODUCTION 


The two-level, fixed effects, full or fractional-factorial design of 
e 2 q)erim.ent, without replication, is the appropriate design for those 
situations where the experiment is very eiqjensive or time consuming. 

An example of costly experimenting is provided by the destructive test- 
mg of simulated aircraft components, as m the rotor burst protection 
testmg described by Mangano (1977). New rotor burst protection mves- 
tigations are planned to measure the containment efficiencies of com- 
posite structures. One such investigation is planned as a two-level 
fractional-factorial experiment. The description (Holms 1977) of that 
experiment is illustrative of one area of applicability of the results of 
the present mvestigation. 

If such an ej^eriment is performed and t observations are obtained 
from t orthogonal experimental conditions, the appropriate emprical 
equation for representing the results can have as many as t terms, 
each with a coefficient that has been fitted to the data. When this is done, 
a question that should be asked is: ”Can the predictive accuracy be un- 
proved if some of the terms are deleted*?" The fact that some of the 
terms might degrade the predictive accuracy of a fitted equation was 
recogmzed by Walls and Weeks (1969) but they gave no procedure for 
identifymg such terms 

A method for the sequential deletion of terms that was mtended to 
reduce the prediction error was given by Kennedy and Bancroft (1971). 
Their method assumed that the e:q)erimenter has a prior established 
order for subjectmg the terms to a sequence of significance tests. Un- 
fortunately, m many experimental situations, there is no basis for estab- 
lishmg a prior order, and m such cases an order statistics procedure is 
appropriate. An order statistics approach for significance testing was 
used m a pair of related papers by Daniel (1959) and by Birnbaum (1959). 
They were not then seekmg to mmimize prediction errors. 

For model selection procedures used with small saturated experi- 
ments (experiments designed to have only as many experimental condi- 
tions as there are model parameters to be fitted), the analysis should 
begm with a mmimum number of estimable terms being sacrificed to 
form a denommator for the test statistic. A procedure using m-terms 
sacrificed, where m can be as small as one, was mvestigated by Holms 
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and Berrettoni (I960) . The procedure was a form, of subset regression 
using backward elimination The object was to delete terms in a manner 
where some control was maintamed over the probabilities of Type 1 or 
Type 2 decision errors. The procedure was called cham pooling and 
used a strategy (iHp, 0 !^, a^) where m^ was the number of terms im- 
tially sacrificed to the test statistic, a was the nominal level of a pre- 

Jr 

limmary test of significance where ’’insignificant” resulted on another 
term (mean square) pooled into the denommator of the test statistic, and 
Q!£ was the nominal level of the fmal test of significance, for the inclu- 
sion of terms in the model . 

The minimizing of prediction error was the object of further inves- 
tigation of the cham poolmg strategy (nip, cUp, or^) as described by Holms 
(1974) . Whereas that mvestigation had assumed that only one cycle of 
analysis would be used, a suggestion given by Holms and Berretom (1969) 
was that more than one cycle of analysis should b© used. The specific 
suggestion was that an analysis be performed with the strategy (mp, ffip, 

Q!j) where mp = 1 to obtain an estimated number of null mean squares, 
fi, tacitly assuming that the population number of null mean squares is 
greater than zero. A second analysis would then be performed with a 
strategy (mp, a^) where mp is an mteger approximation to the 
product r^'^ where r^ had been empirically optimized at about r^ = 0.7. 

The purpose of the present investigation is to use Monte Carlo 
methods to optimize a combined procedure that might contam more than 
one analysis cycle, where the procedure is to be optimized for minimum 
prediction error. Such a procedure would seem to be worthwhile m view 
of the large number of aerospace research and development programs 
where: 

(1) There are a large number of controlled variables 

(2) The experimenting is time consuming or expensive 

(3) The dependent variable has an error variance that might be large 

Based on the results, some specific model selection procedures are 
recommended. 
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PRIOR WORK 


An early practice described by Davies (1956, p. 286) consists of 
pooling some arbitrary number of those mean squares that represent 
higher order terms into an estimate of error variance. When this is 
done, some terms of important predictive value might be deleted from 
the model. Furthermore, such mean squares might unduly inflate the 
estimate of the error variance, thus reducing the sensitivity of the sub- 
sequent procedure, and resultmg m the deletion of additional important 
terms from the model. 

The preservation of sensitivity, when poolmg mean squares mto the 
estimate of error variance, has been an object of the procedure of Daniel 
(1959, andofWilk, Gnanadesikan, and Preeny (1963). 

Daniel (1959) uses the absolute values of the effect estimates as 
order statistics . They are plotted on probability paper and the result 
IS called a half-normal plot. In addition to conditional structuring of the 
AROVA model, Daniel’s objectives mcluded the determmation of "bad 
values, heteroscedasticity, dependence of variance on mean, and some 
types of defective randomization, . . . The half normal plot, com- 
bined with a background of experience, might provide a method by which 
a skillful user could pass judgement on the results of an experiment. 
Daniel concluded that such a plot can be used to make judgements about 
the reality of the largest effects observed only if a small proportion of 
the effect estimates represent real effects. Bimbaum (1959) mvesti- 
gated procedures related to half-normal plotting. His results on . . 
the probabilities of the various possible sorts of errors ..." are 
limited to the single largest order statistic. He surmised that if only a 
small number of effect estimates have nonzero means, the power and 
sensitivity properties will tend to hold approximately. 

The procedure of Wilk, Gnanadesikan, and Freeny (1963) if used 
with 2^ treatments is benefitted if some subjective or prior knowledge 
IS used to decide that tj of the 2^-1 mean squares do not contam real 
effects, or that p = 2^ - tj - 1 mean squares do contain real effects. As 
was shown by Wilk et al. (,1963) , the procedure is not robust agamst 
errors in guessmg the value of rj, and t] must be guessed if the prior 
knowledge is lackmg, which is an assumption of the present investigation. 
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Daniel and Birnbaum have limited their results to experiments where 
only a small proportion of the effects are anticipated to be significant. 

On the other hand, situations can exist where the experimenter might use 
a two-level fractional factorial eigieriment designed such that a large 
proportion of the effects are significant. A need remains for a decision 
procedure that sacrifices only a small number of terms to produce a test 
statistic, and then uses the test statistic to mmimize prediction errors. 

Allen (1971) referred to literature that used the residual sum of 
squares as a criterion for choosing regressors. He pomted out that 
there are at least two objections to using such a criterion 

1. If the residual sum of squares were the only criterion, then all 
of the regressors would be used and there would be no motivation for 
subset regression. 

2 . The residual sum of squares is not directly related to the 
'’natural” loss function which is the mean square error of prediction. 

On the other hand, the methods of Allen obtam a criterion called 
PRESS by predictmg each observation from all the other observations. 
The present investigation is concerned with small saturated experiments 
where each observation is more or less crucial to the estimation proce- 
dure, particularly for predictions at the set of regressor yaiues corre- 
sponding to a given observation. The PRESS criterion is therefore be- 
lieved to be mappropriate to subset regression for small saturated ex- 
periments. 


MULTISTAGE DECISION PROCEDURE 


Assumptions 

The assumptions are as follows. 

(1) The model to be fitted is Imear in the unknown parameters. 

(2) The errors of the observations are mdependently normally dis- 
tributed random variables with a zero mean and a constant variance. 

(3) Orthogonal estimators are available for estimatmg the unknown 
parameters of the linear model. (This orthogonality can be the result of 
the design of the experiment that furnished the observations, or it can 



be the result of an orthogonalizing transformation of the terms of the 
equation, as discussed by Holms (1974) . ) 

(4) An appropriate criterion of the goodness of a subset regression 
procedure is the smallness of the largest of the prediction error mean 
squares over the points of the experiment. 

(5) There is no replication available for an estimate of the "pure" 
error variance. 

The first four of the five preceding assumptions are funadmental to the 
rationale of the method; however, data has been given (Holms and 
Berrettoni (1967)) for believing that the method is robust agamst the nor- 
mality part of the second assumption. The fifth assumption merely ac- 
knowledges the possibility that an altogether different method might be 
preferred m the presence of pure replication. 

The origmai mvestigation of chain poolmg had been concerned with 
three sizes of ei^eriments, namely, experiments furnishing 16, 32, or 
64 observations (Holms and Berrettoni (1967)). The simulations had 
shown drastically reduced decision error probabilities, or the equiva- 
lent, greatly improved information efficiencies, for the larger experi- 
ments. Such results suggest that the method of analysis is relatively 
less critical for the larger experiments, and the methods described 
(Holms and Berrettoni (1967)) are therefore believed to be adequate for 
producmg small prediction errors with experiments providmg 32 or 
more observations. 

For relatively saturated experiments that are smaller than 16 obser- 
vations, the opmion is offered that such eiqperiments are too small to 
provide both (1) good estimates of model coefficients and (2) a good test 
statistic, m cases where random errors are large enough to call for a 
statistical decision procedure. In other words, saturated experiments 
with less than 16 observations should be fitted with models having a fixed 
number of terms with no use of conditional modeling. 

Consistent with the preceding remarks, the simulations of the pres- 
ent mvestigation were all performed with experiments contammg 16 ob- 
servations in the belief that such experiments are large enough to justify 
the use of a statistical decision procedure, but small enough so that the 
precise optimization of the decision procedure would be quite beneficial. 
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Strategy 

The generalization of the strategy (nip, dip, oi^) to be investigated is 
the strategy (nip, Qip 2 , rg, where the symbols are defined 

in appendix A. 

SIMULATION PROCEDUEE 


Population Model 

The experiment is assumed to be a two-level, fixed effects, 
fractional-factorial experiment, where the mdependent variables are 
qualitative or quantitative, and they are controlled to have negligible 
error. In such situations, the mdependent variables are often assigned 
the values of +1 for the ”upper” level and -1 for the 'T.ower'^ level. An 
example of an appropriate model equation for the population mean value 
of the response m the case of four mdependent variables is 

E(Y) = 0!q + + Q!2X2 + ^ “23^2^3 

" 4^4 + ^ 14^4 + “ 24 ^ 2^4 ^IZ4:V2^4: 

“34^3^4 “134^1^3^4 '*'•‘^234^2^3^4 “1234^1^2^3^4 

The subsequent discussion assumes that an equation such as the pre- 
ceding equation will be fitted to the results of a two-level experiment 
where the x’s are ’’design values, " namely, the high level of x^^ is 
represented by Xj^ = +1 and the low level of is represented by 

The smgle observation value of the response is assumed to occur 
accordmg to the model 

Y = E(Y) + e (2) 

where e islhe independently normally distributed random error, with 

E(e) =0; V(e) = (T^ (3) 

The experiment is assumed to be a two-level, fixed effects, 
fractional-factorial experiment where g is the number of mdependent 
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variables, and the experiment is assumed to be a a" fractional repli- 
cate of the full factorial experiment. The observations are assumed to 
result one-for-one from the treatments, and their number is 

t = 2^-^ 


The observations are used to compute estimates of parameters such as 
the parameters of equation ( 1 ) and also to compute mean squares, typi- 
cally by Yates’ method (Davies, 1956, p. 263). 


Assume that the mean squares have been computed in Yates’ 
order and m this order are labeled Z^, Z^, . . . , Z^. The coefficients, 
ck’s, of equation (1) are in Yates’ order. If they are subscripted m that 
order and then given the symbol the expectations of the are 


E(Z^) = 0 -^ + 

1 = 1 , 2 , . . t 


(4) 


O 

The random variables Zjtj are noncentral chi-square variables 
havmg one degree of freedom (Kendall, M. G. and Stuart, A., 1961, p. 
227). Let be the noncentrality parameter. It is related to parame- 
ters already defmed by 



( 5 ) 


A stepwise multiple decision procedure is to be used to delete esti- 
mated coefficients from an equation such as the one illustrated by equa- 
tion (1) . The situation for which the operatmg characteristics of the 
statistical decision procedure will be optimized is assumed to be as 
follows* 

The t single degree of freedom mean squares are assumed to have 
been drawn from t populations. An unknown number, p, of the popula- 
tions (aside from the one associated with /S^) have real effects (X^ > 0 ) 
and the 77 other populations are null populations (rj populations have 
= 0 ). A number, j5, of the populations (other than the one associated 
with jSj^) are to be selected for retention of model coefficients m equa- 
tions like ( 1 ), hoping they have > 0 . 
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The relative magnitudes of the nonzero population parameters to be 
tested can be e:q)ressed by dividing squares of the coefficients by the 
mean value of the'sum of the squares, namely, where the quantities 6^ 
give these relative magnitudes, 


From equations (6) 





k = 1, 2, . . p J 



( 6 ) 


(7) 


Unfavorable Population Model 


The basic cham pooling concept was investigated for the purpose of 
rnirmnizing prediction error as described by Holms (1974), That inves- 
tigation was concerned with the fittmg of models to fractional factorial 
experiments under the condition of population functions of irregular 
shape. The emphasis of the present mvestigation is on the co nditio n 
where the relative values of the population coefficients are all unfavor- 
able to the deletion procedure. For reasons given by Holms and 

Berrettom (1969) this is done by proportionmg the squares of the jS. 

2 ^ 

values to the expected values of the order statistics of a sin^e 
distribution. 


Where the 


are the expectations of the 


order statistics in 

mcreasing order from a sample of size pj the m the order of Yates’ 
computation beyond the first were given a decreasing order by settmg 


^ 2 - 5 - 

°k “ ^p-k+1 
k = 1, . . . , p 


(8) 
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Let X be the arithmetic mean of the noncentrality parameters. 


-Y. 


From equations (5) and (9) 


^ k=l 

2 

Usmg the definition of of equations ^6), 




From equations (10) and (11) 


«2fv 2 

2 


From equations (8) and (12) 






Expectations of order statistics from a gamma distribution with 
scale parameter one, shape parameter 1/2 and many s^inple sizes have 
been tabulated (Harter, H, L. (1964)). Multiplying such values by 2 
gives the expectations of the order statistics of the central distribu- 
tion. Such expectations, for a sample of size p, provide the v^ues 


called for by the definition of the 
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With the distribution of values fixed as just described, the popu- 
lation parameters , as given by equation (14), depend only on the 

^ 2 

single observation error variance c and the mean, X, of the noncentral- 
ity parameters. 

The simulated experiments were performed and the decision pro- 
cedures were mvestigated with t) mean squares havmg A. = 0 and with 
p = 2S-^-l-r, mean squares having the unfavorable set of values 
(eq, (5)) that were just described. 

Observation Simulation 

The number of treatments is 

t = 2^"^ (15) 

In accordance with model equations such as (1) the t observations re- 
sultmg from t treatments are given by equation (2) with E (Yj^) = as 
the variates 


yi = ^l + e^ 

Y 2 = ^2 ®2 

) 




(16) 


In the analysis of an experiment, the application of Yates* method 
estimates each of the t parameters that occur m an equation such as (1). 
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!Furthermore, from equations (12), the nonzero values of the are 


where 


and 



k = 1, 2, 3, . . , p 

I =1, 2 t 

1 = k + 1 for k = 1, . , p (18) 


and l^j = 0 otherwise Let 




X 

gg-h 


(19) 


From equations (17) and (19) 




( 20 ) 
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Any particular strategy (nip, evaluated 

for an array of populations having mp unique values of namely, p^, 
m = 1, 2, . . . , Kip and unique values of the mean noncentrality pa- 
rameter, namely, Xg , f = 1, 2, . . •> 

Because of the freedom to choose the values and because is 

a scale parameter on (eq. (17)) an mvestigation of the effect of varia- 
2 2 

tions in cr IS superfluous, and a will be set equal to one. 

As implied by equation (16) the computer program must perform 
operations equivalent to the following: 

yi = Ml + 6i (21) 

yg = M2 + eg 

^16 = '“16 ^6 (22) 

As previously stated, the are to be computed from the y^ by Yates' 
method. With the bj^ computed, the reversed Yates method of Duckworth 
(1965) IS used to compute predicted values of Y, namely, The predic- 
tion errors are then 

V = - »*l (22) 


Conditions Investigated 


As developed in the subsection on strategy in the mam section on 
MULTISTAGE PROCEDURE, the statistician's strategy consists of the 
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components (m^, r^, ^2’ “f>* 

As developed in^the subsection labeled Unfavorable Population Model, 
nature's strategy consists of the number, p, of non-null population pa- 
rameters, and the mean noncentrality parameter, X. As given by equa- 
tion (19) the relation between the mean noncentrality parameter, X, and 
the scale parameter 0 is 


0 ^=—^ 

2g-h 

For the experunent with 16 treatments, 2^“^ = 16 and the related values 
of 0 and X investigated are as given in table I. 

In general, the smaller the number of null mean squares, rj, the 
greater will be the probability of decision errors. This was illustrated 
m figure 4 of the paper by Holms and Berrettoni (1969). Thus the most 
difficult situation arises for ?] = 0. For the 2®“^ eiqperiment with 
g - h = 4, and where the term (zero degree term) is not subjected to 
testmg, the condition equivalent to = 0 is the condition p = 15. 

A real-life example of a 2^“^ - 2®“^ experiment with p = 15 was 
presented by Daniel (1959) in figure 19 of his paper. As stated by Daniel 
concernmg such a situation: 

"The 15 effects appear to fall mto a nearly normal distribution 
with their own variance, much larger than the error variance. 

The reality of many of these terms is proved by later experi- 
mental work. There may be some sort of central limit theorem 
analog operatmg here on the (discrete) population of effects, but 
no success can be reported in its formulation. " 

One mterpretation of Darnel's discussion is that if an eiqierimenter 
designs his experiment and attempts to select the levels of the controlled 
variables so that measurable effects will be achieved, then a preponder- 
ance of the fitted coefficients of the model equation will have absolute 
magnitudes fairly close to some central value, and a very smaE number 
of coefficients will have much larger and much smaller values, and still 
other coefficients will be scattered between these two extremes, hi 
other words, an approximation to the normal distribution should be ex- 
pected for the coefficients of a model equation resulting from experiments 
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where controlled variables were set with some degree of prior knowledge 
of their effects. 

In the light of the preceding considerations, a nature's strategy with 
p = 15 and a normal distribution of model parameters would seem to be 
a highly likely strategy, and correspondingly, a statistician's strategy 
optimized agamst such a nature’ s strategy should be thought of as an 
approximate Bayes strategy . As developed by Holms and Berrettoni 
(1969), such a normal distribution of model parameters is represented 
by the parameter distributions of table II and these distributions are 
highly unfavorable to the statistical decision procedure. A procedure 
optimized agamst p = 15 and the distribution of of table II may 
therefore also be regarded as an approximate security strategy . Such 
a nature' s strategy (table II) will therefore be chosen as the strategy 
against which the statistician’ s strategy will be optimized. 

Consider the fittmg of a model equation such as that illustrated by 
equation (1) to the results of a two-level, fixed effects, full or fractional 
factorial experiment. The mdependent variables can be standardized so 
that all the x-values are either +1 or -1. With t treatments and t ob- 
servations, the estimators of the coefficients are all of the form 


A 



Vl-" ^ Vt 

t 


(24) 


where the a^, a 2 , . . . , a^ all have values of +1 or -1 The single 
observation error variance of the “Y -values is assumed to be cr^ and 
because the are formed from a linear combination of t indepen- 
dentlY distributed Y-values (eq. (24)), the variance of any is 


V(b ) = t ^ = 2? 
t^ ^ 


(25) 


With equation (1) written as a prediction equation, the predicted 
value of Y, at any one of the pomts of the eiqjeriment, is a linear com- 
bination of t independently distributed b, values, each (eq (25)) with 
variance a /t. With appropriate change of notation, the linear combina- 
tion implied by equation (1) is written 
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( 26 ) 


and from the definition of the x- values m equation 

Ai 

equation (26) are either +1 or -1. Because the 
distributed, and because the are all +1 or ~1, 

V(Y) = V' = 

1 i— , J t 

J=1 


(1), the c^-values in 
are independently 
the variance of is 


(27) 


From equation (27), the reduction of V(Y ) achievable by deleting terms 

Q * 

IS O' A for each term deleted. 

On the other hand, if equation (1) is the population model, and if 
the X- values are all +1, the bias in Y is mcreased by the amount of 
for each value that is deleted. Thus an optimal strategy to mmi- 
mize the squared error of Y should not only delete all terms for which 
the population 13^ is zero, it should also delete at least all terms for 
which the bias contribution to mean square error is less than the vari- 
ance contribution. 


As mdicated by equation (1), for that point of the experiment where 
all of the x-values have the value +1, the expected value of would 



17 


take on its greatest absolute value, which would be 



max 




The values of ^ 0 6^ for p = 15 have been listed m table III. Because 

M j 

cr = 1, these values are also the values of e(y ] jo-. • 

' ^Wax/ 


Reciprocals of EfY,) ja are here defined as coefficients of vari- 

ation for the maximum population mean values. From table III such co- 
efficients vary from a high of 64.3 percent (at 6 =0. 125) to a low of 
4 0 percent (at 6 = 2 000) . This range of such a coefficient of varia- 
tion suggests that the range of 0. 125 gO ^ 2.000 is an adequately wide 
range of 0 to represent the situations that an experimenter might 
encounter 


EVALUATION CRITERIA 

Followmg the selection of terms (where some of the coefficient esti- 
mates are set equal to zero), the predicted values of the dependent vari- 
able can be efficiently computed for all of the combinations of the inde- 
pendent values, by the reversed Yates method of Duckworth (1965). 

Where e are the ’’observation” errors, namely, the pseudo 
om 

normal random numbers generated in the simulation, the "observa- 
tions” are given consistently with equations (16) by 

yoitmn ^ + ®o m (28) 

£ = 1, 2, . . ., m = 1, 2, . . ., m^, i = 1, 2, . . . , t. 

After the model has been fitted and msignificant terms deleted, the 
difference between predicted values, ypugmn’ dependent variable 

for the n^^ simulation and the population mean will be called the predic- 
tion error, and thus it is (consistent with eq. (23)) 
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®pi£inn ^p3jemn " 


(29) 


£ = 1, 2, . . . , £g, m = 1, 2, . . . , m^„ i - 1, 2, , . . , t. 

Over the n simulations, the sample mean square error of predic- 
tion for a given treatment is 


e2 

pi6m „ pifmn 

“e n=l 


(30) 


The maximum of such errors over the treatments is 


4ax = ®^m,max=”=^ (®piim) (81) 

1=1, , , . ,t 

The mean of the squared error over the simulations and over the 
pomts of the space of the experiment is 

t 

®p£m“!^ \_j ®p]£m 
1=1 

Equations (31) and (32) provide two criteria for measuring the effec- 
tiveness of a strategy. The particular set of values of strategy parame- 

—2 

ters that minimizes (as given by eq. ( 31 )) can be called a 

security strategy, and if the pomts of the space of the experiment are 

assumed to be equally likely of bemg of interest, the particular set that 
—2 

minimizes can be called an approximate Bayes strategy. For 

either criterion, the absolute values of squared errors would have been 
the prune consideration. 

An example of a situation where such defmitions of error would be 
appropriate occurs if the expermientor seeks to maximize some pre- 
dicted response, such as the strength of a structure as a function of its 
geometric variables. For such an example, the region of the space of 
the eiqjermient of greatest mterest would be the region m the vicmity of 
the maximum pomt, where the function would most likely have its shaip- 
est curvatures and largest errors due to lack of fit. For such an example, 
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the appropriate criterion to be minimized for the choice of a strategy 

—2 

would seem to be the quantity max equation (31)., 

The criteria of equations (31) and (32) were evaluated using computer 
simulations usmg (m most cases) 1000 experiments. Thus the long run 
mean squared error of the decision procedures was evaluated. This 
leaves open the question of how badly a decision procedure might perform 
m mdividual cases. One approach to this question is to evaluate the sta- 
bility of the mean squared errors observed in the simulations. Thus m 
addition to the criteria of equations (31) and (32) two other criteria for 
the effectiveness of a strategy were mvestigated. They are concerned 
with the stability of the quantities defmed by equations (31) and (32). The 
mstabiLity of these criteria can be measured by the variance of the square 
of the prediction error. If Y is a random variable, the unbiased esti- 
mate of the variance of Y from a sample of size is given by 


r n 


e 


V(Y) = 




n. 






(33) 


The random variate of mterest is the squared error of prediction, namely, 

2 2 
®pjl?mn* equation (33) the estimate of the variance of is 


V 




r n 


e 



\ 


®pifmnj 


n=l 


2 




(34) 


Equation (34) gives an unbiased estimate of the variance of the squared 
error over n simulations. The maximum of this quantity over the 
space of the simulated eiqieriments is defmed by 

(35) 



The arithmetic mean of the variance of the squared error over the 
space of the experiments is defmed by 
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’= 1=1 

Two features not present in the computer program POOLES (Holms 
(1974)) were added to the present program (POOL6U). One of them com- 
putes the average number of terms, selected by the strategy (m^, 

^pl’ ’=’l’ %2’ ■’^2’ each of the values of 0j^, £ = 1, . . •, S-q and 

for each of the population values of m = 1, . . . , m^. The other 
feature computes the ratio of the maximum prediction error to the scale 
parameter 0. The ratio is computed from 9^ and from the e^^^ of 
equation (31) 




_ p£m,max 


COMPUTER PROGRAM 


(37) 


Outline of Program 

Computations were performed accordmg to the FORTRAN-4 program 
POOL6U as listed in appendix B. The antecedents of the program were 
POOL3U (Holms (1966)), POOUVIS (Amlmg and Holms (1973)), and POOLES 
(Holms (1974)). The program is outlined and the parts of POOL6U from 
the earlier programs are identified by the section numbers and titles of 
appendix B in the table that follows. The table shows that only minor por- 
tions of POOL6U were not previously described. The detailed description 


of POOL6U IS mainly limited to sections that are new. 
is given m appendix C. 

Illustrative output 

Section number 

Section title 

Prior program 

lA 

DECLARATIONS AND TABLES 

POOLMS 

IB 

INPUTS AND CONSTANTS 

POOLES 

1C 

POPULATION MEANS 

(new) 

ID 

STRATEGY 

(new) 

2 

SIMULATIONS AND MODEL 

POOL3U 


FITTING 



Section number 
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3 

4 

5 

6 

7 

8 
9 


Section title 


Prior program 


CONSTRUCTION AND ORDERING POOL3U 

OF MEAN SQUARES 

DELETION OF TEEMS (new) 

PREDICTIONS POOLES 

ACCUMULATION OF ERRORS POOLES 

DETERMINATION OF MAXIMUM POOLES 

AND MEAN SQUARED ERRORS 
OUTPUT (new) 

YATES' METHOD SUBROUTINE POOLES 


Details of Program 

Section lA, - Declarations and tables . - The mput of the critical 
point values (the statistical tables) is the same as in POOLMS (Amling 
and Holms (1973)). The declarations are similar to POOLES (Holms 
(1974)). The NAMELIST output was incorporated to facilitate checking 
of the program and was not otherwise used. 

Section IB. ~ Inputs and constants . - The constants defining the 
populations, the experiments, and the sequential deletion strategy, are 
read from data cards m the following order, with the order of the fields 
bemg the same as the order of the s 3 nnbols m the following description. 


Format Description 

(13A6, A2) REMARK (I), arbitrary literal information such as 

particular use of program, date of last change, and 
so forth. 


(415) LGMH, NE, MRHO, KODE 


(I8.6F8.3) LTH, (THETA(L), L=l, LTH) 

14/ KTEM 

(10F8 . 5) (DELTA(K, M) , K=l, KTEM) 


22 


Format Description 

There are as many (10F8.5) cards as are necessary to 
read (DELTA(K, M),K=lj KTEM). Furthermore, the 
card sequence for KTEM and DELTA (K,M) is repeated 
as many times as necessary to satisfy the statement 
(appendix B) "IK) 5 M=l, MRHO." 

(415, 2F5 . 3) MP, KP(1) , KP(2), KF, RETA{1), RETA(2) 

(The associated HEAD statement is actually m 
Section ID.) 

Section 1C. - Population means . - After the initial constants have 
been read, the next major operation is the formation of the population 
mean values . The number of population sets to be exammed durmg the 
mvestigation of a strategy is the number resultmg from all combmations 
of the number, m , of p-values, and the number, ofD-values. 

fj {/ 

With respect to equation (1) all the population model parameters are 
first set equal to zero with the DO-ioop endmg at Statement 10. The non- 
zero values of the B are mitially set equal to 5 using the DO-loop 

-1* -1-^ .ill 

endmg at Statement 20. The DO-loop endmg at statement 20 serves the 
purpose of equation (20) with (7=1 and 0=1. The value of cr = 1 is 
retained, but the adjustment for 0 is made after the population mean 
values have been computed. 

With the population jS- values (aside from 0) established at state- 
ment 20, the objective is to compute the population mean values from the 
/3-values by the reversed Yates' method (Duckworth (1965)) . The first 
step is to reverse the order of the /S-values, which is completed at state- 
ment 22. The use of the reversed Yates' method then yields the array 
YMUM (I,M) as completed at statement 40 for all the values of p, 

YMUM (I,M) IS therefore an array of population means This 

array of population means is to be ejq)anded over the mean noncentrality 
parameters, , to give the effect of equation (20) . This effect is pro- 
duced on the population mean values by the multiplication 
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and this operation is completed with the creation of the array TMU (I, L,M) 
at statement number 49. 

The values of YMU are thus indexed over treatments i, i = 1, . . . , t 
over arbitrary values of 0^; S. =1, . . . , : and over the different values 

of m = 1, . . . , m^. 

Section ID. - Strategy . - In terms of mathematical symbols previously 
defmed, the strategy parameters are functions of numbers that are read at 
Statement 50 as follows: 

Argument Function 

FORTBAN symbol mathematical symbol 


MP 

m 

P 

KP(1), 

“pi 

RETA(l) 


KP(2) 

“p2 

RETA(2) 

^2 

KF 

“f 


For each strategy investigated, the contents of the error arrays are 
first set equal to zero with.the DO-loop ending at statement 99 . 

Section 2. -• Simulations and model fitting . - This section is essen- 
tially the same as that described for POOLES by Holms (1974), except 
that m POOL6U the experiments are simulated and fitted as if they are 
full-factorial experiments. 

Section 3. - Construction and ordermg of mean squares . - This sec- 
tion IS essentially the same as that described for POOLES by Holms 
(1974). 

Section 4. - Deletion of terms . - This section contams the major 
changes between POOL61J and POOLES. The major distinction is the 
provision for more than one iteration of sequential deletion, the maximum 
number of iterations being three. The new strategy parameters are r^^. 
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rgj Qf^g, The flow chart for the iteratioas is shown by figure 1. 

The statements from ”DO 415 J=1,MPX" to ”420 JETA = J-1” are 
essentially the same as those m POOLMS (Amlmg and Holms (1973)) 
from "DO 15 J=l, M" to ”20 JETA = JA-1. ” 

_ The statements beyond 420, and up to and mcluding statement 424 
provide for the revision of the number of mean squares imtially pooled 
m the given sequential deletion strategy in accordance with the value as- 
signed to r^ or Tg. 

The statements beyond 424 and termmatmg with statement 429 set 
the msignificant model parameters equal to zero, as was done m state- 
ments 416 through 419 of POOLES (Holms (1974)). 

Section 5. - Predictions . - This section is an abridgment of that 
described for POOLES (Holms (1974)). 

Section 6. - Accumulation of errors . - This section is the same as 
that of POOLES (Holms (1974)). 

Section 7. - Determmation of mayrmmn and mean squared errors . - 
This section is the same as that of POOLES (Holms (1974)) except for the 
computation of the arithmetic mean, p, of the number of teinns (beyond 
the zero degree term) that are retamed by the decision procedure, and 
the computation of COEBMX. 

Section 8. - Output . - The output is illustrated in appendix C . The 
NAMELIST output was incorporated only for program checkmg. 

Section 9. - Yates* method subroutme . - This subroutme is essen- 
tially that of part of the mam program of POOLMS (Amlmg and Holms 
(1973)) except with the last few statements modified so that the subroutme 
can be used for the direct Yates’ method and also for the reversed Yates’ 
method; as was also done m POOLES (Holms (1974)). 

The program POOL6U contams some statements from POOLES that 
are not essential to the present purposes of POOL6U. 
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SIMULATION RESULTS 


Monte Carlo Sample Size 

Resists of a preliminary investigation of the effect of Monte Carlo 
sample size on the stability of the empirical results are shown by figure 2, 
hi general, the results converge to a constant when the number of sampled 
experiments is 1000 or more. Variability of results occurs as the number 
of sampled ei^eriments is reduced below 1000. All of the strategy com- 
parisons were performed for 1000 sampled experiments. 


Large Coefficient of Variation 


One of the conclusions of Holms (1974) was that a widely useful three 
parameter strategy of sequential deletion for mmimizmg prediction errors 
consists of the strategy (m^, a^, = (1, 0.75, 0.10). The possibility 

of performmg two analyses in sequence had been suggested by Holms and 
Berrettoni (1969). That suggestion had included performing an initial 
analysis with the strategy (in.p, 0 ;^, with m^ = 1. The values of 0 ;^ 
and would be selected accordmg to some attempt to control the proba- 
bilities of type 1 or type 2 errors. The mitial analysis would yield an esti- 
mate of Tj, Then based on table 3 of Holms and Berrettoni (1969), a second 
analysis would be performed with m^ approximately equal to 0.67 or 0. 75 
times the rj from the first analysis. The joint implications of the preced- 
mg two sets of results are that for the present mvestigation with six param- 


eters, a useful strategy might be as follows, (m , a r,, a 

P px X P2 




“f ) = 


(1, 0.75, 0.700, 0.75, 0.700, 0.10). Some Monte Carlo experunentmg was 
started with strategy parameters m the immediate vicmity of the set just 
listed. 


The most adverse value of p among nature’ s strategies is the value 
p = 15 Using this value of p, the simulations were performed for all 
the listed values of 6 of table I. 

The strategies were most often compared in terms of the maximum 
coefficient of error, defined at equation (37), The largest values 

of ^aa TYi occurred at 6=1. 125 . The better strategies for 0=0. 125 
are listed m table IV which also lists the values of for both 

0 = 0.125 and 0 = 2,000. 
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The best strategies of table IV for C„. (0.125) for each value 

60 ^ JJLl^ 

of m are listed m table V together with additional information on the 
P 

operatmg characteristics of these strategies . These strategies are all 
good strategies at d = 0.125, which represents the situation where the 
model parameters are relatively small (table III) in comparison with the 
error standard deviation (a = 1. 0). (From previous discussion in the 
section on simulations and from table III, the ratio of the maximum popu- 
lation value of the dependent variable to the standard deviation at 0 = 
0.125 IS 1.556 whereas that ratio at 0 = 2.000 is 24.898.) 


Small Coefficient of Variation 

If the statistician's loss function is the maximum relative error 
over the space of his experiment, then the strategy that is 

optimal for 0=0, 125 is a security strategy because withm the present 
mvestigation 0 is larger for 0 = 0.125 than for any other value 

(table I) of 0 mvestigated. 

On the other hand, if the statistician’ s loss function is simply the 

absolute value of the maximum squared error over the space of his 

—2 

experiment, namely, e„„„ as defmedby equation (31) then that quantity 
IS a ma xim um withm the present mvestigation at 0 = 2. 000 and thus the 
security strategy for such a loss function would be the strategy that mim- 
mizes Cgg (2.000). The strategy that mmmiizes ^ee max 
IS the strategy with no deletion which is symbolized as (mp, r^^^, 

“p2* ^t] 2’ » » ’ Some other strategies 

that were good with m^ = 1, 0 ;^^^ = 1. 00, and ~ 1. 00 are listed with 
resultmg values of C.. (0.125) and (2.000) in table VI. 

003 XXL/v 003 111. A 


Admissible Strategies 


For the purposes of the present mvestigation, a strategy will be 
classed either as admissible or as dommated according to its values of 
„^(0) at both 0 =0.125 and 0 =2.000. A strategy will be said to 
be dominated of for 0=0. 125 there is another strategy with the same or 
lesser Cgg (0 . 125) and with a lesser mx 
will also be said to be dommated if there is another strategy with the 


same or lesser C 


ee,mx 


(0.125). 


(2.000) and with a lesser C 
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Any strategy that is not dominated is defined as bemg admissible Por 
example, witnm table V, the strategies for = 1, mp = 3, and mp = 5 
are admissible, whereas the other three strategies are dominated 

In addition to the strategies Listed m tables IV and VI for 0 !p^ = 1.00 
and Q!£=:1.00, many strategies with smaller values of and 
were mvestigated. Among aU of the strategies mvestigated, the admissi- 
ble strategies, together with some of their operating characteristics are 
listed m table VII. 

Security Regret Strategy 

The strategies of table VII mclude the strategy mp = 0, which re- 
sulted in the smallest observed value of (2.0), namely, 0.5282, 

and they mclude the strategy (5, 1.00, , 0.05, 0.675, 1.00) which 

gave the smallest observed value of C (0.125), namely, 7.689. 

The regret function of a statistical decision procedure , as a function 
of a parameter 6 , is here defined as the excess loss occurring with the 
procedure at a particular value of d as compared with the loss that 
would have occurred had the best statistical decision procedure been used 
for that particular value of 6 . For the purposes of the present investiga- 
tion a regret function E(0) is defmed for 6 = 0.125 as being the 
(0. 125) for any strategy divided by the value of for the best 

' ' illJC 

strategy for that value of 0, namely, 7.689, and R(0 ) is defined for 
6 = 2 . 000 as being the C (2. 000) for any strategy divided by the 

value of for the best strategy for that value of 0, namely, 

00 y IIlX 

0.5282. 

Thus, 


R(0.125) =Ceg^^^(0.125)/7.689 

and 


R(2.000) = (2.000)/0.5282 

00 y XllX 

The smgle strategy that has the smallest regret function over both 
0=0. 125 and 0=2.0 is defmed as the security regret strategy . The 
security regret strategy is thus the sequential deletion procedure which 
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produces the least increase in prediction error for p = 15 and an un- 
favorable distribution of parameters over that prediction error which 

could have been achieved if the best strategy had been chosen for the 

2 

given (unloiown) value of error variance, cr . 

Examination of the R(0) values of table VII shows that the parame- 
ters of the security regret strategy (mp, Op^, r^^, a^ 2 > ^2’ “f^ 

1.00, , 0.75, 0.360, 1,00) for which the values of E.(0) are 1.0966 

for 0 =0.125 and 1.0682 for 0 = 2.000. 

Selection of a Strategy 

hi summary, if the experimenter wishes to minimize the maximum 
prediction error over the points of the ei^ermient when the variance 
error is relatively large (0=0. 125) the statistician’ s strategy should be 

(5, 1.00, , 0.05, 0,675, 1.00), which minimizes C (0.125), 

(tables V and VII). If the experimenter wishes to mmimize the maximum 
prediction error over the pomts of the experiment when the variance 
error is relatively small (0 =2.000), the satistician should use the 
strategy mp = 0, which is to say that no model deletion will be performed, 
which minimizes (2.000) (table VII). 

If the eiqperimenter has no basis for a choice between the two pre- 
ceding extreme strategies, a security regret strategy with (mp, Qipj^, r^^, 

Q!p 2 , ^ 2 , ttf) = (1, 1.00, , 0.75, 0.360, 1,00) should be used, in 

which case the values of the regret fvmction (from table VII) will be E.(0) = 
1.0966 if 0 =0.125 and R(0)= 1.0682 if 6 = 2.000. These values of the 
regret function show that the relative prediction error standard deviation 
will be mcreased by at most 9.7 percent over what it would have been if 
the worst value of 0 had occurred and the best strategy agamst it had 
been used. Thus the security regret strategy of (mp, oip^, r^, o-pg, r 2 » 

ctj) = (1, 1.00, , 0.75, 0.360, 1.00) must be concluded to be a widely 

useful strategy. 

Variance of Predicted Squared Error 

The three strategy selections described in the preceding section are 
based on a Monte Carlo mvestigation that reported mean values of predic- 
tion errors over 1000 simulations. The quoted results thus teH what the 
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mean long run results will be as a function of strategy selection. The 

subject of short run results was not discussed. Some insight into the 

short run performance can be gamed by examining the observed values 
2 

of V(e This quantity gives the observed variance, for samples of 

size 1000, of the maximum squared prediction errors over the simula- 
tions, as defmedby equation (35), IE this variance is relatively small, 
then operatmg characteristics such as are relatively constant 

XXX 

2 

from simulation to simulation. But if V(e )_„ is relatively large, then 
the short run performance of a strategy could be erratic. 

In the case of large coefficients of variation (small values of 6) the 

2 

strategy performance was not erratic - the values of V(e were 
small for all the strategies of table VII for 0=0. 125, The strategy per- 
formance can be erratic for small coefficients of variation (large values 

0 

of 6). Thus the values of V(e were large or small for 0 =2.000, 
depending on the strategy (table VII) . This response to 0 shows that the 
bias component is the component of the prediction error that can be 
erratic . 

2 

The response of V(e ) to the strategy parameters is to be exam- 
med, and this will be done for 0 =2. 000. First of all, the parameters 
and r^ are relatively unimportant. Namely, exammation of ta- 
ble VII shows that identical operatmg characteristics are often obtained 
despite wide variations of and Furthermore, the three specific 

strategy recommendations of the preceding section all use = 1. 00. 

P“** 2 

For these two reasons, the influence of strategy parameters on V(e ) 

' it* ^ 

was exammed only for cases with 0 !^^ constant at oip^ = 1. 00. Thus the 
only strategies of table VII that were examined for their influence on 

9 

V(e )^^ were those that mitiated the first analysis cycle as specified by 
TTip and 02^2 then did or did not perform a second analysis cycle as 
specified by r£ and 

Usmg results from table VII with 0 =2, 000 and = 1. 00, an arbi- 
trary second degree polynomial was fitted. The form of the polynomial 


was; 
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In 


= <*0 ^'l“p + + V 2 + ^4“f + ^“p 

* ^6“p2 + Vi + ^8“f + ^9“p“p2 + '^lomprj + 


* %2“p“f 


* %sV“t * ®14’'2“f 


(38) 


The fitting used the backward deletion procedure of Sidik (1972) -with a 
nominal significance level of 0.01. Results rounded to two significant 
figures are as follows; 

In [v(e^)jJ = 9. 1 - 0.44 mp + 5.5 0!p2 

- 6.7 Qfj - 11.2 Q!p2 + 9-1 r2 + 1.0 ^p“p2 


The coefficients of equation (39) show that the variability of 

Cp.p(2.000) ^y decreases with increasing mp, decreases rapidly with 

increasing 0 !p 2 J increases rapidly with r 2 , increases with the product 

of nip and a ^2 decreases with increasing a^. The range of nip 

included the values m = 1, 2, 3, 4, 5, Thus even for m„ = 5 the 

p p 

cross product term m mpO!p 2 together with the first degree term in 

Q!p 2 9-1*6 of slightly less influence than the second degree term m Q^p 2 - 


The values of the strategy parameters all serve to control the num- 
ber of terms deleted, which fixes the number, p, of terms retained. 

Thus a given value of p can result from many different combinations of 
values of the strategy parameters. From the results exhibited by equa- 
tion (39) the conclusions as to what combination of strategy values would 
result in a given p while mmimizmg V(e^)^^ are. mp is unimportant, 
ap 2 should be small; r 2 should be small; and should be large. 

Thus to minimize both and V{e )^J a strategy should 

use = 1.00 and the smallest values of Q'p 2 and rg that give accept- 
able values for over appropriate values of 0. The choice of 

06 iJQJK. 

m IS not critical, and the use of a , = 1. 00 is generally acceptable, 

P ^ 
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CONCLUSIONS 

An mvestigation was conducted to determine what statistical tech- 
niques should be used for model fitting to the results of a two-level, 
fixed-effects full- or fractional-factorial, orthogonal experiment with 
16 treatments. Multiple sequential deletion strategies involving as many 
as two preliminary tests and one fmal test were evaluated, using Monte 
Carlo techmques, under the criterion of mmimum prediction error. 

Three strategies were identified as bemg appropriate depending on 
the extent of the experimenter’s prior knowledge. 

1. If the experimenter has prior knowledge that the relative error is 

relatively small (coefficients of variation of 4 percent or less) then no 
deletion should be used. The strategy is (m^, 0 :^^, r^, ~ 

> ) 3 ) > » ) • 

2. If the experimenter has prior knowledge that the relative error 
is quite large (coefficients of variation m the neighborhood of 65 percent) 
the strategy should immediately delete the five smallest absolute value 
terms and then test with contmued poolmg at a nominal test level of 0. 05 
to estimate a number rj of msignificant terms. The number of terms 
deleted from the model should then be the integer value of 1 + 0. 675 77 . 

The strategy is (m^, r^, o-^^) = (5, 1.00, 0.05, 

0.675, 1.00). 

3. Without prior knowledge about coefficients of variation, a security 
regret strategy should be used. It consists of beginning the sequential 
deletion with the single smallest absolute value term sacrificed to the 
test statistic. Testing is done with continued poolmg at a nommal test 
level of 0.75 to estimate a number tj of insignificant terms. The num- 
ber of terms’ deleted from the model should then be the integer value of 

1 + 0.36 The strategy is (m , a r,, a 2 > cUf) = (1, 100, , 

0.75, 0.36, 1.00). 
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APPEM)IX A 
SYMBOLS 


Mathematical rORTRAN 

symbol name 


Description 


B(I) 

Cee,mx COEKMX 

E (. . .) 

e RN(I) 


estimate of j3^ 

ratio of maximum prediction error to 
scale parameter, eq. (37) 

e:q)ectation of . . . 

single observation random error 


max 


ERSQMX 


S 

g-h LGMH 

h 


i,j,k, I,J,K 

Ji,m, n L,M, N 

KF 

KODE 

KP 


£ 


9 


m 

P 


LTH 


maximum over treatments of mean 
square prediction error over simu- 
lations 

number of mdependent variables 

experiment contams 2^“^ treatments 

experiment contams (1/2) times num- 
ber of treatments m full factorial 
experiment 

subscripts 

mdex number for «£, Amlmg and 
Holms (1973) 

amount of NAMELIST output 
desired 

mdex number for 0 !^, Amlmg and 
Holms (1973) 

number of 6 values mvestigated in 
any computer run 

number of mean squares pooled before 
testmg begins 


MP 
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Mathematical 

FORTEAN 

Description 

symbol 

name 



MRHO 

number of p values mvestigated m 
any computer run 

% 

NE 

number of simulated experiments in 
any computer run 


RETA(l) 

number of mean squares pooled for 
testmg at level a ^2 mteger 

value of 1 + VjT], fig. 1 

^2 

EETA(2) 

number of mean squares pooled for 
testmg at level is mteger value 

of 1 + rg^, fig. 1 

t 

IT 

number of treatments 

V(. . .) 


variance of . . . 

V(e^) 

^ 'max 

7ESQMX 

maximum over treatments of sample 
variance of mean square prediction 
error over simulations, eq. (35) 



mdependent variable 

■^1 


conceptual value of dependent variable 


YOBS(I) 

observed value of dependent variable 

2i 

Z(I) 

mean squares m Yates’ order 

“f 


nominal significance level of fmal test 

“pi 


nommal significance level of first pre- 
liminary test 

“p2 


nommal significance level of second 
prelimmary test 


B(I) 

regression coefficients m Yates' order 


DELTA(K,M) 

parameter determmmg relative magni- 
tudes of coefficients m population 
model, eqs. (6) and (8) 
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Mathematical, 

symbol 

rOKTEAN 

name 

Description 



ejqjectation of order statistic of 

2 

of a variable 

1) 


number of mean squares having non- 
centrality parameter of zero 

A 

V 

ETA 

number of mean squares concluded to 
be null during any analysis 

H 

THETA(L) 

scale parameter 

X 


mean over eiqjermient of noncentrality 
parameters, eq. (9) 

\ 


noncentrality parameter 

i\ 

YMU(I, L, M) 

fVi 

population mean value of for i 

treatment 

p 

KTEM,KRHO(M) 

number of non-null coefficients {num- 
ber of values) beyond zero order 

term m population model 

A. 

p 


number of coefficients concluded to be 
non-null on any analysis iteration 

■p 

AVRHO 

mean number of coefficients concluded 
to be non-null m a strategy mvestiga- 
tion 

<j 


standard deviation of e 



n n n 
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APPENDIX B 
PEOGRAM LISTING 


c 

c 


TECLARATIOMS AND TABLES 

DIMENSION REHARKU'Of ALPHAtll), TB{6H,1C>, RNtlfe), IWO« 1 6 ) , Z ( 1 6 ) ♦ 
1TPETA(6), KRH0(3J» DELTA(16,3Jt YMuH{16,3)t YHU(16t6i3l, 
ZAVRh 016,3), ER50(l6t6 ,3), EPSOS 0 1 16 , 6 , 3 I , ERS CMX 1 6 1 3 ) , COERHX « 6 , 3 J , 
3AVERSQ(6i3) t VES0HX(6,3), AVV£SO(6,3>t KPCZli PETA(?), ALFA(2> 

CCHMON KK, Y0eStl6)t ®C16) 

DATA CALPHA(H,I=l,ll)/Ci.Oei,C.OC!2,''.005,r.01,0,L25,O.D5,0.10f 
IC.ZS^f.SO.O.TS,! ,r/ 

DATA T(TB(I,JJ,J = l,l£.)tI = li2‘»>/O.D,D.OiD.'',O.C,C.C,C.C,O.P,0.afC.O 
1,C.D,2.00D0,1.999<>9,1 . 9<»99 7 , 1 . 999S6 , 1 .99 9 17 , 1 .99 687,1.98 77,1.923.1 
2.7C6,1.382»?.9976,2.9 96L,2.99CM,2.9Pu9,2.951,2.9«.«,2.8C6,?,527,2 .; 
38 6, 1.688, 3.976,3.962, 3 .®29 , 3 . P7o, 3 . 76 P, 3 .625, 3. M i? ,2 . 9«|9 , 2 . 395 , 1 .9 
9 6 1,4.887, 4, 8 1»5, A. 758, 9.65, 9. 99, 9. 21, 3. 89 , 3. 28 7, 2. 658 , 2. 189, 5. 79, 5. 

563. 5. 96. 5. 31. 9 . 99. 9 . 68, 9. 28, 3. 57, 2. 893,2 . 37 1,6. 51, 6. 33, 6. 11, 5. 67, 5 
6. 96, 5. G9, 9, 61, 3. 83, 3. 1 1,2. 59, 7. 22, 6 . 96, 6 .65, 6 .35, 5 . 88, 5. 99, 4.9 1 ,9. 
70 6, 3. 29, 2. 6 0,7. 81. 7. 5 2, 7. 1C, 6. 78, 6. 26, 5. 79, 5. 17, 9. 27, 3. 45, 2. 82, 8. 3 

89. 8, Cl, 7. 53, 7. 17, 6. 59, 6. 03, '.4 1,9. 95, 3. 6L, 2. 95, 8. 82, 8. 99, 7. 95, 7. 53 

9 .6.89.6.28 .5.61 .9 .62, 3.79 ,3.67,9.26 ,8.89 ,8.33 ,7.87 ,7. 13,6 .5C.5 .8 1 , 
A9 . 77, 3. 87, 3. 17, 9. 67 , 9. 21, a. 68, 8. 16, 7. 37, 6. 71, 5.99, 9. 92, 3. 99, 3. 27,1 
BC,D5,9.Ss,8.95,P.92,7.59,6.91,6.15,5.CS,9.1D,3.37,ir.9P,9.e6,9.2r, 
C« . 6 6, 7. 7 9, 7. 07, 6. 30, 5. 17, 9 . 20, 3. 96, 1 :.72,1C. 19, 9. 93, 8. 83, 7. 96, 7. 23 
0 ,6.99,5.29,9 .3u,3.5', 1 1.01 ,lo.4P,9.69,9, CC,8. 12, 7. 38 ,6.57,5.9C ,9 .3 
£9,3.63,11 .20,10.64,9. 89,9.17,8.28,7.52,6 .69,5 .5L,9 .9 a, 3.7", 11 .S’ ,1 
FC.8 6,10.u3,9.39,8.93,7.65,t.61,5.6r,4.S6,3.77,j.l.76,ll.r>,lC.22,9. 
651, 8. 58, 7. 78, 6. 92, 5. 69, 9. 69, 3. 89, 11. 98, 11. 28, 1C. 9C, 9. 67, 8. 72, 7. ?r, 
H7.C 3, 5. 76, 9. 71, 3. 90, 12. 19, 11. 98, 10. 55,9.83, 6. 86, 8, 22, 7. 13, 5. 87, 9. 7 
IP, 3. 96, 12. 39, 11. 65,10.76, 9. 99, 8. 5 9, 4. 13, 7. 2 3, 5. 95, 9. as, 9. '’2, 12.se, 
Jll.e7,lC.93,lC.19,9.12,fi,29,7.33,6.r3»9.92,9.:s/ 

DATA (t TB II, J), J= !, 10), 1 = 25, 971/12, 76, 12, 11. IS, 10. 29, 9. 23, 8. 39, 

17 , 9 2, 6, 11, 4. 98, 9. 19, 12. 93, 12. 22, 11. 26, ID. 93, 9. 39, 8. 99, 7. 51, 6. IP, 5. 
209,9.19,13.09,12.38,1 1.91,1 r. 5 6, 9. 4 9, 8. 5 9, 7. 6 ",6. 2 5, 5. 1C, 9, 29, 13.2 

39. 12. C3, 11. 55, 1C. 68, 9. 59, a. 63, 7. 68, 6. 32, 5. 16, 9. 30, 13. 39, 12. 68, II. 6 
98, 1 C. 78, 9. 64, 8. 72, 7. 76, 6, 38, 5. 22, 4. 35, 13. 53, ir. 82, 11. 8C, 1C. 88, 9. 79 

5. 8. 81 , 7. 83, 6. 99 , 5. 2P, 9. 90, 13.67, 12. 96, 11 . 91, 1C. 98,9.83, P.89, 7. 9r, 6 
6. PC, 5 . 33, 9, 95, 13. 80, 13. C9, 12. Cl, 11. '>7, 9. 91, 8. 97, 7. 97, 6. 56, 5, 38, 9. 5 
70, 1 3. 93, 13. 21, 12. 10, 1 1 . 16, 9. 99, 9. 09, 8. C9 , 6. 62 ,5.93 , 9 , 59, 19 . C 5, 1^.3 

82. 12. 19, 11. 25, 10. 17, 9 . 11, 8. 11, 6. 6 8, 5. 9 8, 9. 5 8, 14.17,13 .43,12.27,11. 
934, 10. 15, 9. 18, 8. 17, 6. 79, 5. 53, 9. 62, 19. 29 ,13. S3 ,12. 35,11. 93, 10. 22,9. 
A 2 5, 8. 23, 6. 8", 5. 58, 9. 66, 19. 4 1,1 3. 6 3, 12. 9 3, 11. 5 1,1 0.29, 9. 3 1,8. 29, 6. 8 
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I 

R5,5 .63,4 .7G, 14.53« 13.73«12*51,IU59 ,10.36|9»37»8 .35*6.90,5.67,4.74 
C, 14 .64, 13. 8?, 12.59 , U .67, 10»43, 9.43, S.4J ,6.95 ,5.71,4 .78,14.75, 13.9 
01 ,12. 67,1 1.75, 10. 5 0,9. 49, 8. 46 ,6. 99,5.75, 4 .82,14. 85,14. 00, 12.75, 11. 
E8 3, 1C. 57, 9. 55, 8 .51,7.03,5 .79 ,4 .86 , 1 4 • 95 , 14 .09 ,12.83,11 .90,10.64,9. 
F61, 8. 56, 7. 07, 5. 83, 4. 9 0,15. 05, 14. 17, 12. 90, 11. 9 7, 10. 70, 9. 67, 8. 61, 7.1 
Gl,5 .8 7, 4, 94, 15. 15, 14. 25, 12, 9 7, 12. 04, 10. 76, 9, 72, 8. 66, 7, 15, 5. 91, 4, 98 
H, 15 .24,14.33,13.05, 12 .11, ID, 82, 9. 77, 8. 7 1,7. 19 ,5. 95, 5. 01,15.33, 14-.4 
I 0,1 3. 12, 12. 18,1 0.88, 9,82,8, 76, 7. 23, 5. 99, 5, 04, 15.42,14.47, 13.19,12. 
J2 5, 1C. 94, 9. 87, 8.81 ,7.27,6.03,5.07/ 

OAT A ( (TB 1P> , 1=48,641/15.50,14.54,13 .26,12.32 ,11.00,9.92, 
18 .8 5,7.31 ,6.07,5.10,15.58,14.60,13,32,12 .38,11.06,9.97,8.89,7.35,6 
2. 11 , 5. 13, 15. 66, 14. 66, 13. 38, 12*44, 11. 11, 1C. 02, 8. 93,7. 39, 6. 14, 5 . 16,1 
35.73,14.72,13.44,12.50,11 .16,10.07,8.97,7.43,6.17,5.19,15.80,14 ,79 
4 , 13 .5r,lZ,56, 11 .21 ,10.12,9.01,7.4 7,6.20,5.22, 15*87,14.85,13.56,12. 
5 6?, 11. 26, ID, 17, 9. 05, 7. 5 1,6. 23, 5. 25, 15. 93, 14. 9 1,1 3. 62, 12, 6 8,11. 3 1,1 
6 0.2 1,9.09,7.55, 6.26,5 .28, 15.99, 14.9 7,13.67,12 .73,1 1.36,10.25,9, 13, 
77 , 5 9, 6. 29, 5, 31, 16. 05, 15. 03, 13. 72, 12. 78, 11.40, 10. 29, 9. 17, 7. 63, 6, 32, 
8 5 ,3 4,16.11,15.10,13.77,12.83,11 .4 4, 10.33,9,21 ,7.67,6.3 5,5.37, 16,17 
9, 15 ,16,13.82,12 .88 ,11 .48,10.37,9.25,7.70,6.38 ,5.40,16.23,15.22,13. 
A87, 12. 93, 11. 52, 10, 4 1,9. 29, 7. 73, 6. 41, 5. 43, 16, 2 9, 15. 28, 13. 92, 12. 97,1 
Bl. 5 6, 10, 45, 9. 33, 7. 76, 6. 44, 5. 46, 16. 34, 15. 34, 13. 97, 13. 01, 11 .60,10.4 9 

C, 9. 37,7.79,6.47,5.48,16.39,15.40, 14 ,02,13.05,11.64,10.53,9.41 ,7.82 

D, 6. 5C, 5. 50, 16. 4 4, 15. 46, 14. 06, 13, 09, IX. 67 ,10 .5 7, 9. 4 5, 7, 85, 6. 5 3, 5. 52 
F,C.C,r.C,0.0,C.C>,0.D,C.C,0.D,D. 0,0. 0,0.0/ 

C 

NAHELIST /OUTl/ M /0UT2/ YOBS /0UT3/ YHUH /0UT4/ H,L 

NAMELIST /0UT5/ YMU /0UT6/ N /0UT7/ RN 

NAMELIST /0UT8/ J.INDX 

namelist /0UT9/ MPX,KPX,KF 

namelist /OUT 10/ B 

NAMELIST /OUTll/ IND,Z,e 

NAMELIST /0UT13/ N,M,L /0UT14/ J,JN, TEM /0UT15/ ETA 

NAMELIST /0UT16/ NC 

NAMELIST / /OUT20/ ETA.AVRHO 

NAMELIST / /0UT22/ ERSO ,FRSOSO 

C 

C IP.- INPUTS AND CONSTANTS 

C 

READ(S,8rO) (REMARK II ) ,1 = 1,14 ) 

WP1TE(6,801) <REMARK(I),I=1,14) 
rEAD(5,8021 LGMH, NE, MRHO, KOOE 
KG= LRMH 
1T= 2**LGHH 

WRnE(6,8D3) KG, IT, HRHO, NE 

ITM1= lT-1 

ITM2 = IT -2 

FIT= TT 

FITHl = ITMl 

FNE= NE 

IF (NE.EO.O) FNE = l.D 
NEMl = NE - 1 
FNE HI = NEMl 

IF (NFKl .EO. 0) FNEMl = l.D 
REA0(5,sn4» LTH, (THETA(L1, L=1 ,lTH» 

DO 5 M=1,HRH0 

READI5,806> KTE H , ( DEL T A (K ,M 1 , K=1,kTEM) 

KRHO(M»= KTEH 
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WR1TE(6,807> KTEM^ I DELTil I K » W ) , Krl.KTEH) 

S CONTIK^UE 
IPr ?«^ckg 
FIP r IP 
IPpl= IP + 1 
FENPIP r FNf^FIP 
FFK lir r F^JFMli^FTP 
C 

1C,- POPULATION KEANS 
KKz 

nr ^0 Mri,MpHo 

z ^ 

DO 1C IrljIP 
Bp) z C 
1C CPNTIfUE 

KTEHr KRHO {«) 

KTEhP;^ - KTFM + 2 
DC 21 I-1,KTEK 
F Cl + nr DELTA (I^M) 

2*> CCNTI^UE 

UODE tGT, D ) WRITE <6, OUT 1 ) 

22 T::i,IP 
if'PlMI- IPPl - T 
YCBS(1>= P(IPt^lMl) 

^2 CONTINUE 

(K'^DE .GT, I ) WRITE (6, 0UT2 ) 

CALL YATES 
DC t: I- 1 , ip 

IPPIKT- IPPl - I 
Y”U^MT,P):: P(IPP1KI1 

continue 

4T CONTINUF 

IF (KODE .GT. 2 ) WRITE (6t OUT 3 ) 

DO 49 M-1,PRH0 
!• DX F z H 
DC 4 8 L = 1 , L TH 
INDXL z L 

I*" (KCDE .GT. 3 ) WRITE (6, 0UT4 } 

Dr 47 1 = 1, IP 

YPUn,L,^*> = YHUH (I,H )»THET Aa) 

47 CCNTJrUE 
4H CONTINUE 
49 CONTINUE 

IF (KODE .GT. 4 > WRITE (6, OUTS ) 

10. - STRATEGY 

READ <5, 608, END = 899 ) HP , K P ( 1 ) ,KP ( ? > , K F ,RET A C U , RE T A ( 2 ) 
DC 55 NC=1,2 
KPX= KP(NC) 

ALFAr#C)= AtPHACKPX) 

55 CCNTJNUF 

CALL Sand (xs) 
nr 99 M=1,HRH0 
DC 98 L=1,LTK 
AVRHO(L,K)= C*D 
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rc 97 lrl,lp 
ERSO(I,L,M>r 0.? 

EPSCSO* I ,L,K) = D.D 

97 CGNTIMUE 

98 Cf’^^I^UE 

99 rCSTINUF 

2. - sthulattons and hcofl fitting 

I'P 699 n=i,NE 
I^DX^’ = N 

ir (KODE .GT, 5 > WRITE C6, 0UT6 ) 

IF (Fr.E.0.0) GO TO 2;;C 
GC 213 T-1«IT 
TALL '’ANOCRNd)) 

215 CONTINUE 

IF (K'^DE .GT, 6 J WRITE (6, 0UT7 ) 
r.O 21*^ 1 = 1, IT, 2 
F= SGRT{-2.F*AL06CRN( 1 >>» 
r= 6.2831853*RN CI + 1 ) 

pr ( I )= E + Cns tD) 

RKd + i):: E*SlN(nj 

21G CCNTINUF 

IF «KCDF .GT, 6 > WRITE (6, 0UT7 ) 

6C TO 201 
2’~r Pr, 229 1 = 1 , IT 
RN (I J = •'.0 

Z'-9 CONTINUE 
Zri DO 69' H=il,NRHO 
INDXN = “ 

DC 6SC L=1,LTH 
ITDXL = L 
KK= LPMH 

IF (KODE ,GT. 12) WRITE (6, 0UT13) 

2 ’'^ pc 21A 1 = 1, IT 

YPSS«1)= YMU(I,L,M) +PN(I) 

214 CONTINUE 

IF (KCDE .GT. 1 ) WRITE (6, OUT? ) 

CALL YATES 
GO TO 3DD 

3, - CONSTRUCTION ANT ORDERING OF MEAN SQUARES 
C 

3'v OP 3:9 1=1, IT 

INO(D= I 

Z (I>= B(I+11*B(I+1 )/FlT 
Bin = B(T) / FIT 
3CR CONTINUE 

IF (KODE .GT. 10) WRITE (6, OUTll) 

IF (NP ,LT. 1) GO TO 432 
C 

DC 313 J=1,ITM2 
TEST= ZdTMl ) 

IN= ITNl 

00 312 NA=J,ITH2 
IF(TEST-Z(NA)) 312,312,311 
3ll TEST = Z(NA> 
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IN= NI 

312 CONTINUE 
ITEHz INOCINI 
TEHr ?(IN» 

IKDtlNir IND(J) 

ZdNiz Z(JI 
IKD(J1= ITEM 
Z(J} - TEH 

313 CONTINUE 

IF (KOOF .6T, ID) WRITE <6, OuTll) 
4 DELETION OF TERHS 
MPX= WP 

IF CKODE .GT. 12) WRITE (6, OUT13) 
DC 424 NC = 1,2 
IKDXNC = NC 

IF (KCDE ,6T, 15) WRITE (6, OUT16) 
KPX z KPINC) 

IF (KGDE .GT, 8 ) WRITE (6, 0UT9 ) 

IF (KPX ,EQ. 11) GO TO 424 

MFPlz HPX + 1 

JKz MPPl 

TE«= C.O 

DC 415 J = 1,HPX 

INDXJ z J 

TFh = TEH + Z(J) 

IF (KODE .GT. 13) WRITE (6, 0UT14) 
4 15 CONTINUE 

no 419 J=rHPPl,ITHl 
INDXJ z J 
FJN r JN 

TEST - F JN*Z( J)/(TEM ♦ 2CJ)) 

IF (TFST - TP(JN,KPX)> 416,416,420 
416 TEH = TEH + Z(J) 

JK z JN + 1 

IF (KODE .GT, 13) WRITE (6, 0UT14) 
419 CONTINUE 

JFTA z iThI 
GC TO 422 
42C JET‘A = J-1 
422 ETA = JETA 

IF CKCDE .GT. 14) WRITE (6, 0UT15) 
IF (KODE .6T. 12) WRITE (6, 0UT13) 
IF IRFTA(NC) .LT. 1.0) GO TO 421 
MPX z JETA 
GC TO 424 

421 MPX z 1 ♦ IFIX(RETAINC)*ETA) 

424 CONTINUE 

IF (KODE .GT, 8 ) WRITE (6, 0UT9 I 

MPPl : HPX ♦ 1 

JN z KPPl 

TEH z 0.0 

DO 42E J=l, HPx 

INDXJ z J 

INDX z IND(J)+1 

B(INDX}z 0.0 
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7EK = TEH + Z(J) 

IF <K0DE .6T. 7 I WRITE (6, 0UT8 ) 

IF <KODE .GT. 9 > WRITE (6, 0UT1D> 

4 25 CONTIK'UE 

IF ( KF .EO. 11 » GO TO 434 
DO 429 J=HPP1,ITM1 
IROXJ = J 
FJN ~ JH 

TEST = F JN*Z { J) /ITEM + 2(JI) 

IF (TEST - TBCJNjKF)) 428,428,430 

428 I^DX = INDCJ) + 1 

IF (KODE .GT. 7 1 WRITE ( 6 , OUTB ) 
P(INDX) - O.D 

IF (KODE .GT. 9 ) WRITE (6, OUTIC) 

429 continue 
ETA = FTTHl 
GO TO 433 

4^j JFTA = J-1 
ETA r JETA 

IF ( KODE .IT. 8 » GO TO 433 
IF (KODE ,GT, 14) WRITE (6, 0UT15) 

DC 431 J=JETA,1T 

INDXJ c J 

INDX = T\D(J>+1 

IF (KODE .6T. 7 ) WRITE (6, 0UT8 1 

4^1 CONTINUE 

IF (KODE .GT. 9 I WRITE (6, OUTIO) 

GO TO 433 

432 ETA = n.O 
GO TO 433 

4 34 E7A = HPX 

433 AVRHO(L,H) = AVRHO(L,H» ♦ FITHl - ETA 

5. - PPEDICTIONS 

5 CO KK - KG 

540 DC 546 I-1,TP 

IPPIMI r IPPi-I 
YOBS(T) - BdPPlHI J 
546 CONTINUE 

IF (KODE .GT. 1 1 WRITE (6, OUT? > 
CALL YATES 

IF (KODE .GT. 9 ) WRITE 16, OUTIO) 

GO TO 600 

6. - ACCUMULATION OF ERRORS 

6CC DC 619 Irl,IP 
IPPIHT - IPPl-I 

TEH= (B(IPP1HI> - YMU(I,L,M))**2 
Epson, L,H) = ERSO{T,L,M) ♦ TEM 
EFSCSS(I,L ,MJ= FRSOSQ ( I,L,H) ♦ TEM**2 
6C9 CONTINUE 
6BD CONTINUE 
690 CONTINUE 

IF (KODE .GT, 14) WRITE (6, 0UT20) 

IF (KODE .GT. 21) WRITE (6, 0UT22) 
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IF (MF .EO. 0» GO TO 7C0 

6 99 CONTINUE 

?♦- DETE9NINATJ0N OF MaXIHUH ANQ HEAN SQUARED ERRORS 

TRD DC 79F M=1,HRhO 
DC 7b: L=1*LTH 
C = C.D 
D = C.D 
E= C.r- 
F= C.C 

DO 75C 1=1, IP 
C = APAX1(C,ERSC(I,L,M>1 
0 = D + ERSO(I,L,M) 

TEM= ERSOSQ(I,L,M)- { I ERSQ (1 ,L , M > ) **2 ) /F NE 
E = AMAXI (E,TEH) 

F r F + TEH 

7 50 CONTINUE 

EPSOMX<L,HI = C/FNE 

CO£RHX(L,H> = (SQRT<ERSOHX(L,M) ))/THETA(L) 
AUERSO(L,H» = D>/FENBIP 
IF <NE ,E0. 01 GO TO 78P 
VFSOHX«L,M)r E/FNEHl 
AVVESO<L,H)= F/FEHITP 
AURH0{L,HT = AVRHO ( L , H 1 /FNE 
782 CONTINUE 
7^0 CONTINUE 


8.“ OUTPUT 


WRITE (6,809) 
WRITE (6,811) 
WRITE (6,813) 
WRITE (6, SI'S) 
WRITE (6,817J 
WRITE (6,831) 
WRITE (6,817) 
WRITE (6,83S) 
WRITE (6,817) 
WRITE (6,836) 
WRITE (6,817) 
WRITE (6,837) 
WRITE (6,817) 
WRITE (6,838) 
WRITE (6,817) 
GO TO 5C 
8®9 STOP 
C 


(MP, ALFA(l), RETA(l), ALFA(2), RETA(2), ALpHA(KF)) 
(KRHO (M ) , H=1 ,HRH0» 

(THETA(L),(AVRH0 (L,M» ,H=1 ,HRH0 >,l = i,lthi 
( THETA! L ) ,(ERSQMX (L, H » , M z 1 ,FRhO ) ,L= 1 ,LTHJ 
{ THFTAtL ) ,(AVERSO (L , H » , H= 1 ,HRHO ) ,L=l ,LTHI 
( THE TA( L) ,(VESOHX(L ,H) ,M = 1 ,MRH0 ),L = 1 ,LTHI 
( THETA( L ) , (AVVESO(L,M) ,H=1 ,HRHO ),L = 1 ,LTH) 

( THETA(L),(COFRHX (L ,M) ,H=1 ,MRHO >,L = 1 ,LTH) 


SCO FORMAT (13A6,A2) 

8C1 FORMAT ( IH 1 , // lOX , 1 3A 6 , A2/ / ) 

8C2 FORMAT (9T5) 

803 FORMAT (1H0,3X,4HKG =I5,5X,4HIT =15 ,5X .6HHRH0 =I5,5X,4HNE =15) 

804 FORMAT (I8,6F8.3) 

8F6 FORMAT ( 1 4/ ( 1 OF 8 . 5 ) ) 

807 FORMAT «1H0,5HRH0 =15 , 5X , 7HPEL T A =// (IX , lOFlG .5 ) ) 

808 FORMAT (415, 2F5.3) 

8C9 FORMAT ( IHl // , 1 X , 3HMP r 1 5 , 5 X , 9H ALPHAP 1 =F 6 .3 ,5X,7HR ET A 1 =F6,3,5X, 
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S 11 
813 
8 15 
817 
6-^1 
8'^5 
8 36 
837 
8 38 


19HflLPWAP2 rfb.3,5X ,7HRETA2 
FORPAt (1H0,5MRH0 =3113//! 
FORMAT ( 1H0,5HTHETA) 

FORMAT ClH0,2CXr‘>HAVRHO//) 
FORMAT { 1X,F8.3»3E14.4) 
FORMAT UH0,20X,6HERS0MX//) 
FORMAT ( lH0t20X,6HAVERS0//) 
FORMAT ( 1H0,20X,6HVFS0MX//) 
FORMAT (lHc»2CX,6HAVVESQ//) 
FORMAT ( lH0,2DX,6HC0ERMX//> 


=F6.3,5X»8HALPHAF 


END 


=F6.3) 


SUBROUTINE YATES 

R.- YATES METHOD SUBROUTINE 

COMMON MK,Y I16J »B (161 
II r 2*#KK 
IIDB2 = II/2 
KKMl = KK-1 
DO 908 K=1,KKM1 
DO 906 li:l»II,2 
IP102 - (l-»lI/2 
B(1P1D2) r Y(l-*l)-*Y(Il 
LL = IP1024IIDB2 

906 B(LL) = Y(I+1)-Y{IJ 
DO 907 1=1,11 

907 Y {II = B(1 I 
90B CONTINUE 

DO 909 1=1,11 ,2 
1P1D2 = (I+11/2 
B(1P1D2! = Yd^D+Vm 
LL = IP1D2-»I1DB2 
BfLLI = ¥tI+ll-Y(I» 

909 CONTINUE 
RETURN 
END 




AUG 2A, 1977 IF 

INEHl WAS IF (FNEHI 

IF mtiA 

SUBSCpIPTEO JULY 29, 

1977 


K6 = 

A IT = 16 

HUHO = 3 NE - 

1000 




fTHO * 11 

DELTA = 






1 •98710 

1.519S0 1.2A600 

1.04300 .S7S50 

•72960 

.59850 .47710 

•3626C 

•25260 

RHO = 13 

DELTA s 






2«C*>2&0 

•30S30 

1.59640 1.33540 

.21630 .12370 

1.14090 .96200 

•84470 

.72200 .60950 

•5C44Q 

• 404 60 

RHO - 15 

DELTA = 






2«10619 

•43665 

1.66452 l.4C9'<9 

.3S211 .26965 

1.22186 1.06945 

.16921 .10835 

•93855 

.82213 .71606 

•61764 

.52503 


illustrative output 



HP- I ALPMA^I 1 .000 RfTAl z .000 AIPHAP2 

RHO - 11 13 15 


THETA 

AVRHG 


.12 5 

.1566+02 

.1371+^2 

.1371+02 

.2sr 

.1 368+02 

.1370+02 

.1372+02 

.5CC 

.1370+C2 

.1375 + 02 

.1367+02 

1 ,COC 

.1379+02 

.1375+02 

.1368+02 

'2.C0C 

.1381+02 

.1 385 + 02 

.1372+02 


ERSOfiX 


.12 5 

,1112+ni 

.1112 + 01 

.1111+01 

.25 0 

.1105+01 

.1106 + 01 

.1110+01 

.5DC 

.1100+01 

.1116+01 

.1112+01 

l.OOC 

.1099+01 

.1096+01 

.1120+01 

2 . DO 0 

. 1116 +ri 

.1112+01 

.1273+01 


AVERSO 


.125 

.1005+01 

,1005 + 01 

.1005+01 

.25C 

,1006+01 

,1007 + 01 ' 

.1007+01 

.5CC 

.1008+01 

.1010 + 01 

.1018+01 

1 .OOC 

.1006+01 

.1010+01 

,1095+01 

2.00C 

.1009+01 

.1011 + 01 

,1 129+01 


VESOKX 


.125 

.2309+01 

,2284 + 01 

.2276+01 

.25 D 

.2290+01 

.2318 + 01 

.2301+01 

.5HC 

.2302+01 

.2317+01 

.2401+01 

1 .coc 

.2289+01 

.2310+01 

.2513+01 

2.CQC 

,2306+01 

.2335 + 01 

.382^+01 


RFTA7 


360 


alphaf 


1. ClO 








^ivveso 


.125 

• 2*>C •198a-^01 

.50C .2GD6+D1 

l*OOC •2nCjO+Cl 

2.C'*»c .i<?en4^i 

COERf*X 


.12t .e'43'l+Ol 

•25C •42U5+C1 

.501 .20se'^01 

1. coc .io<4e + (i 

2. C0C .52fi2+ro 


.19ft6+01 .1992+01 

.1992+^1 •1997+01 

.2015+01 .2r55+01 

.2004 + 01 .2165 + 01 

.2024+01 .25B7+C1 


.8438 + Pl .8432+f/I 

.4237 + 01 .4214 + 01 

.2113 + 01 .21C9 + D1 

.1047 + 01 .IrSfi + Ol 

.S272+no .5642+PC 




cn 
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TABLE I 

INVESTIGATED VALUES OE SCALE PARAMETER 
AND ASSOCIATED VALUES OF MEAN 
NONCENTRAUTY PARAMETER 

Scale parameter, Meaa noncentrality parameter, 


e 

X 

0.125 

0.25 

.250 

1.00 

,500 

4.00 

1,000 

16.00 

2,000 

64.00 
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TABLE II 


UNEAVORABLE DISTRIBUTIONS OF RELATIVE 
VALUES OF MODEL COEFFICIENTS 
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TABLE m 


PARAMETEB COMBINATIONS, ^yr=0dy^, p-- 15 


0.125 


0,250 j 0.500 


1.000 


2.000 


2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 


2.1082 

1. 6645 
1.4094 
1. 2219 
1.0694 

.9386 

.8221 

.7161 

.6176 

.5250 

,4368 

.3521 

.2699 

.1892 

.1084 


0,2635 

.2081 

,1762 

.1527 

.1337 

.1173 

.1028 

.0895 

.0772 

.0656 

.0546 

.0440 

.0337 

.0237 

.0136 


0,5270 

.4161 

.3524 

.3055 

.2674 

.2346 

.2205 

.1790 

.1544 

.1312 

.1092 

.0880 

.0675 


1, 0541 

.8322 

.7047 

.6110 

.5347 

.4693 
.4110 
, 3580 
.3088 
.2625 

.2184 

.1760 

.1350 


.0473 I .0946 

I 

.0271 i ,0542 


2.1082 

1. 6645 
1.4094 
1.2219 
1.0694 

.9386 

.8221 

.7161 

.6176 

,5250 

,4368 

.3521 

.2699 

,1892 

.1084 


4.2164 

3,3290 

2.8188 

2.4438 

2.1388 

1.8772 
1. 6442 
1.4322 
1.2352 
1.0500 

.8736 

.7042 

.5398 

.3784 

.2168 
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TABLE V 

ADDITIONAL PROPEBTIES FOB STRATEGIES 
THAT WEEE GOOD AT 0 = 0. 125 

[Q!p^ = 1.00, Q!£=1.00, p = 15, iig = lOOO.J 



































TABLE VII 

PEOPERTIES OF STRATEGIES ADMISSIBLE FOB 0=0 125 AND FOR 0 = 2 000 


[B(0 125) = 689, R(2 000) = C^^^^O 5282, p ^ 15, - 1000 ] 
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TABLE VII - Continued 
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FigurB 1 - Strategy (iHp. r^, a^2/ '^2< “f> 
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(c) EFFECT OF Hg ON Cgg 

Figure 2 - Effect of sample s{ze on Monte Carlo results 
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